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SUMMARY 


The local recovery factor was determined exjierimentally along the 
surface of a thin-walled 20° included angle cone for Mach numbers near 
6.0 at stagnation temperatures between 1200° R and 2600° R. In addition, 
a similar cone configuration was tested at Mach numibers near 4.5 at 
stagnation temperatures of approximately 6l2° R. The local Reynolds 
number based on flow properties at the edge of the boundary layer ranged 
between O.lXlO^ n.riB 3.5X10^ for tests at ten^jeratures above 1200° R and 
between 6x10^ and 25X10^ for tests at temperatures near 6l2° R. 

The results indicated, generally, that the recove 2 y factor can be 
predicted satisfactorily using the square root of the Prandtl number. 

No conclusion co\ild be made as to the necessity of evaluating the 
Prandtl number at a reference teng)erature given by an engilrical equation, 
as opposed to evaluating the Prandtl nunber at the wall temperat-ure or 
static temperature of the gas at the cone surface. 

For the tests at ten 5 )eratures above 1200° R (indicated herein as 
the tests conducted in the slip-flow region) , two definite trends in the 
recovery data were observed - one of Increasing recovery factor with 
decreasing stagnation pressure, which was associated with slip-flow 
effects one of decreasing recovery factor with increasing temperature. 
The true cause of the latter trend could not be ascertained, but it was 
shown that this trend was not appreciably altered by the sources of 
error of the magnitude considered herein. 

The real-gas equations of state were used to determine accurately 
the local stream properties at the outer edge of the boundary layer of 
the cone. Included in the report, therefore, is a general solution for 
the conical flow of a real gas using the Beattle-Brldgeman equation of 
state. The largest effect of temperature was seen to be in the terms 
which were dependent upon the intemal energy of the gas. The pressure 
anB hence the pressure drag terms were unaffected. 
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INTRODUCTION 


The transport of energy across the hoimdary layer siir rounding bodies 
flight has heen the subject of intensive theoretical and experimental 
st\idy. A significant number of the studies have been concerned with the 
amount of available energy absorbed by the bcdy; hence the recovery- 
factor term has become considerably important. Although various theories 
have been developed for the prediction of the recovery factor, a lack 
of experimental data exists for high Mach nunbers at the correct flight 
enthalpy levels. 

The various theories as developed can be divided into two main A 

groups. The first of these is applicable at low speeds (Mach numbers 
below 3 • 5 )* In these theories (see refs. 1, 2 and 3) air was considered 
to be a perfect gas with thermodynamic and transport properties invariant 
with temperature. This assumption led to the results of Crocco (ref. l) 
who found that the recovery factor for a laminar boundary layer was 
predicted by the square root of the Prandtl number. A conclusion of this 
analysis was the Invariance of the recovery factor with flight Mach 
number and Reynolds number. 

The second approach was necessitated by the advent of flight speeds 
greater than Mach numbers of 3.5 wherein the assumption of constant 
thermodynamic and transport properties was no longer valid. Several 
authors (refs. 4 through Y) , upon examining the basic energy transport 
processes, found, first, that the enthalpy rather than the tenperature 
of the gas must by used in evaluating the recovery factor; secondly, the 
recovery factor for a laminar boundary layer could still be predicted by 
the square root of the Prandtl number, provided the Prandtl number is 
evaluated at a proper reference temperature. This approach then accounted 
for the real-gas effects (up to dissociation) and indicated that the 
recovery factor would be a function of both temperature and Mach number 
(see ref. 4 ). 

The present investigation provides data showing the local recovery 
factor along a 20° cone at a free-stream Mach number near 6.0 and stag- 
nation temperatures up to about 9^ percent of flight stagnation teiipera- 
tures. The data obtained are compared with tne square root of the Prandtl 
number. The Reynolds number range of the tesns conducted at high 
teirperatures falls near the middle of the region defined by Tsien and 
others (see ref. 8) as the slip-flow region; consequently, the relations 
which can affect the recovery factor in slip ’low are used. In the 
continuimi region some data for a 20° cone at 'fech numbers near 4.50 and 
temperatures of about 6 l 2 ° R are presented. 

In the analysis of the high-temperat\ire data, it was believed 
necessary to have accurate Information conceiving the conical flow field. 
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Consequently, a general solution for the conical flow using the Beattie 
Bridgeman equation of state was developed to provide the needed 
information, and this solution is included herein (see appendix A). 


NOTATION 


A 

a 

Cf 

Cy 

h 

K 

Kn 

k 

k* 

I 

M 

m 

Pr 

P 

Q 

Re 


area, ft^ 

speed of sound, ft/ sec 
skin-friction coefficient 
dimensionless heat-transfer coefficient 
specific heat at constant pressure, ^ Btu/llD °R 
specific heat at constant volinne^ Btu/lh ^ 
enthalpy Btu/lb 

correction term (defined in eq. (B23) or in eq. (AlO) 
Knudsen number 

thermal conductivity^^ Btu/ft hr 
dependent variable given in equation (9) 
length of cone , ft 
Mach number 

model vail thickness perpendicular to model axis^ ft 

Prandtl number, — ^ 
k 

pressure^ 2 lb/ft2 

heat transfer rate Btu/hr 

Reynolds mmiber^ 


^Without subscript m^ the values for Cp^ and p pertain to 

air. 

Without subscript t;, the values for p^ h, and p denote 
static conditions 



R 


gas constant^ ft-lb/lb 


r 

T 

t 

U 

X 

X 

y 

7 

e 

Q 

b 

P 

T 

9 


c 

m 

n 

s 


recovery factor, -- ^ (for Cp asstimed constant) 

absolute temperature^^ ^R 

temperature 

velocity^ ft/sec 

mole fraction of particular species in a mixture of gases 
distance along model longitudinal axis^ ft 

coordinate perpendicular to x axis^ ft 

c 

ratio of specific heats ^ 

Gy 

emissivlty 

radius vector angle from axis of symmetry^ deg 
characteristic temperature of vibrations^ ^R 
viscosity^ Ib/ft-hr 

dependent variable given in equation i! 9) 

density^ 

time^ hr 

cone half angle ^ deg 


Subscripts 

conditions at outer edge of the boimdfry layer along the cone 
surface 

model material 
normal to radius vector 
shock vave 


^“^^See footnotes 1 and 2 p* 3* 
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t total conditions (i*e.^ conditions that vould exist if the gas 

were hroui^ht to rest isentropically) 

tr transverse cuivature 

V conditions at the vail or surface of the model 

w free-stream conditions 

Q a particular radius vector 

( ) conditions employing zero flov and nearly zero pressure 

ZiT 

Superscripts 

^ slip flov; ideal gas flov in appendix A 

' quantities evaluated at the reference temperature (see eq. (6)) 

DESCRIPTION OF APPARATUS 
Wind Tunnels 


In this investigation tvo different vlnd tunnels vere used to test 
tvo similar models. For the tests conducted at high temperatures^ herein 
referred to as the hot-flov tests ^ a pehhle-bed heater vas employed in 
conjunction vith the Ames Lov-Denslty Wind Tunnel. Where the stream 
stagnation tenp>e natures vere just high enough to prevent liquefaction 
of the air components the tests vere made in Ames 10- by l4-inch tunnel 
and are referred to as the cold-flov tests. 

Hot-flov tests . These tests vere conducted in a Mach number 6.0^ 
lov-density^ high-temperature vlnd tionnel (fig. l) vhich had as a heat 
source zirconium-oxide pebbles (dla. = 3 / 8 - inch) heated by natural gas 
burners. The pebbles vere alloved to reach a maximum temperature of 
kOOO^ F. Surroimding the pebbles vas an inner liner of insulating 
refractory vhich vas contained in a steel shell fabricated in several 
sections. At the bottom and top of the shelly openings vere provided 
for the entrance of the cold-storage air. A circular steel nozzle vas 
placed in the center section of the steel shell. The nozzle^ vhose 
Internal contours vere obtained by the method of characteristics for a 
cold-flov Mach number of 6.0 (boundary- layer corrections applied)^ had 
an exit diameter of A. 296 Inches and a throat diameter of 0.468 inch; 
and vas I 5 . 2 O inches long. A vater jacket surrounded the entire nozzle 
to keep the inner vails cool. The model vas placed in an open- jet-type 
test section. Dovnstream of the test section a converging -type diffuser 
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was used to decelerate the flow. The air, upon leaving the diffuser, 
was cooled by a water injection system, and the desired hack pressure 
was regulated hy a five-stage steam-ejector system. 

Cold-flow tests .- The Ames 10- hy li<-Inch Supersonic Wind Tunnel 
was used for these tests. A complete description of the tunnel is given 
in reference 9* 


Models and Instrumentation 


Hot-flow tests .- The model used in these tests was a 20° included 
angle cone fabricated from type 4l6 stainless steel. A sketch of the 
model is shown in figure 2(a) . The configuration was mounted on a 
translte end-pli;ig locked to a hollow stainless- steel sting. (The 
transite material was used for the end-plug to insure negligible heat 
conduction losses out the base of the model.) Care was taken to prevent 
nonuniform circumferential heat conduction by maintaining a uniform wall 
thickness (the wall thickness was 0.025 in.). In order that a reliable 
value of emlsslvlty could be used over the temperature range encountered, 
the external and internal surfaces of the model were plated with a thin 
coating of pure nickel and then "flashed" with platinum (ends of 
thermocouple were covered by the plating) . 

Chromel-alumel thermocouples were embedded along the surface of the 
model and in the radiation shields which are shown in figure 2(a). 
Chromel-alumel was used because of its proven perfoimance for measuring 
ten5)eratures over a range from -300^ to 2500^ F. The thermocouple leads 
(vrire size vas No* 28 B and S gage) vere ‘brought out the sting and 
connected to indicating microvolt potentiometers. To insure consistent 
readings throughout the investigation^ the model thermocouples were 
calibrated before and after each test against known controlled teirpera- 
tures. During the tests the temperature -time history of the cone surface 
teii5)eratures was obtained by photographing the indicating dials on the 
potentiometers. A 70-mm motion-picture camera capable of 20 frames per 
second was used with an electric timer placed in view of the camera. 

A traversing mechanism inside the test cJiamber provided means for 
momting several models and/or probes and mov .ng these individimlly into 
the center of the air stream. The model for these tests was always 
placed at zero angle of attack with respect to the geometric center line 
of the nozzle. 

The total temperature of the hot-flow air stream was measured by an 
aspirated double- shielded chromel-alumel thermocouple probe shown in 
figure 3- The thermocouple leads from the probe were connected to a 
standard potentiometer which was initially balanced against a standard 
cell. Concerning the measurement of the stagnation ten5)erature^ it was 
realized that effects such as radiation^ conduction^ and boimdary- layer 
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phenomena along the probe could cause significant errors in obtaining a 
true stagnation temperature reading. The probe was designed to prevent 
these effects from influencing the true stagnation readings (see ref. 10 
for a discussion of these types of temperature effects) ^ and the 
thermocouple was calibrated against known temperatures. 

The free-stream static -pressure in the test section (assumed constant 
across the test section) was measured at two stations. One station^ a 
wall orifice^ was located l/2-inch upstream of the nozzle exit. At the 
other station, two probes (alined perpendicular to the nozzle center 
line) were located just aft of the nozzle exit and out of the nozzle 
flow. The pressures at the wall orifice and one exit probe were recorded 
on Pirani gages. The other exit pressure was measured with a McLeod 
gage (see ref. 11 for a description of these gages). The stagnation 
pressure of the entering flow was measured by a standard U type mercury 
manometer. 

Cold- flow tests .- The model used in the 10- by l4-inch wind tunnel 
was also a 20^ included angle cone fabricated from type ^l6 stainless 
steel (fig. 2(b) ) . The model was sting-mounted in the tunnel test 
section at 0^ geometric angle of attack. 

The several teii 5 )eratures along the surface were also measured with 
chrome 1-alumel thermocouples embedded at the various stations along the 
model (see fig. 2(b)). The thermocouple leads were connected to indicating 
microvolt potentiometers and the system was calibrated in the same 
manner as was indicated in the previous description of hot-flow tests. 
Instrumentation for determining total temperature and Mach number has 
been discussed in reference 9* 


TEST PROCEDURE 
Hot -Flow Tests 


The tunnel was operated between temperatures of 1200^ and 2600^ R 
at stagnation pressures of approximately 6.6 and 12.4 psia. For each 
test condition the flow was assumed stabilized when the static-pressure 
data recorded by the one McLeod and two Pirani gages were identical. 

The Mach number was determined from the ratio of static- to total- 
pressure data. (Corrections were applied for the effect of caloric 
imperfections on the ratio of static pressure to total pressure.) Also, 
the Mach number was checked during several high-teirperature runs by a 
flow visualization method in which the angle of the oblique shock 
originating at the cone apex could be obtained within ±0.2 . For all 
runs the average Mach number was found to be 5*92 with a maximum 
deviation of ±0,12. 
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When the flow was established, the shielded thermocouple (total- 
temperature probe) was placed in the center of the air stream, and the 
total temperature was recorded. The probe waij removed and the model 
was rotated into position. The temperature o;;' the model was then 
measvired continuously until it remained constant with time at all stations 
( in lire© surface and invo radiainion shield inec^erainures) • Upon reaching 
equilibrium (approximately 10-13 minutes) the model was removed from the 
hot stream and the total ten 5 )erature probe again was placed in the hot 
stream. The temperature of the hot stream was taken as the average of 
the two temperature probe readings. The difference between these 
readings was never greater than 60 F (a function of temperature 
distribution in the heater) , but for 80 percent of the runs the 
difference was less than 30 ° P. 

To determine the cone external radiation losses to the surroundings, 
the model was allowed to cool under conditiont. of zero flow and nearly 
zero pressure. For these tests the model was brought up to equilibrium 
temperature. Then, the tunnel pressure was quickly reduced to a pressure 
of less than 10 microns of mercury and the flew was shut off. Data were 
recorded continuously up to the time at which the model was within 50 ° 
to 100° F of its environmental temperature (usually a time Interval of 
approximately 20 minutes) . It was observed during these tests that the 
internal radiation was approximately zero since the temperature of the 
outer radiation shield was approximately the same as the temperature of 
the cone (see appendix B for the derivation of the radiation heat loss). 


Cold-Flow Tests 


In the 10- by 1^-inch wind tunnel the data were obtained at Mach 
numbers between 4.4l and 4.59 (in conjunction with a different investi- 
gation) and the total temperature of the air, heated electrically, was 
of the order of 612 ° R. The Mach number was calibrated and adjusted as 
indicated in reference 9* It was believed that these data would be 
representative of cold— flow data obtained at M — 6.0 in the continuum 
flow region since the difference in Mach nimiber should have no significant 
effect on the recovery factor. 

In these tests, after several minutes the model reached equilibrium 
(assimied zero external and internal radiation Losses) and the potentiom- 
eters vhich indicated the surface temperatures vere photographed. 


REDUCTION OP DATA 
Hot -Flow Tests 


Gas pr operties .- In the analysis of the at.nodynamlc heat-transfer 
measurements made during the present investigation, it was necessary to 
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A 


determine the local stream properties at the outer edge of the boundary 
layer of the cone. Because of the high ten^jeratures encountered during 
the present tests^ the gas can no longer be considered thermally and 
calorlcally perfect (ideal) and real-gas equations of state must be 
used; thus a general solution for the conical flov of a Beattie-Bridgeman 
gas was developed* This solution is Included in appendix A, and for 
conparison puiposes table II show's the ideal and realigns properties 
(up to dissociation) at the outer edge of the boimdary layer for several 
cone angles at various tenperature levels* 


- ? factor .- As mentioned previously, the significant quantity 

eflnlng the recovery factor for high-temperature flows is the enthalpy 
of the air rather than the temperature. Thus the recovery factor is 


r 


hp-h 


c 



( 1 ) 


The enthalpy terms h^, and h-t^ refer to characteristics outside 

the boimdary layer and in the free stream, respectively, ami were 
calculated from the flow solutions developed in appendix A and from the 
wind-tunnel flow properties. The quantity h^ is the recovery enthalpy 
and is defined as the enthalpy of the air at the surface under the 
conditions of no heat transfer to or from the surface. The latter 
condition did not exist during the present Investigation. Furthermore, 
the skin temperature of the cone rather than the enthalpy of the surface 
air was measured. It was necessary, therefore, to develop an involved 
set of equations to reduce the skin-temperature data to the desired 
recovery enthalpy. The development and discussion of the equations is 
presented in appendix B. 

Briefly, the equation for the recovery enthalpy developed in 
appendix B can be written 


hp == Iv + K (2) 

where K may be considered a correction term, and accounts for the 
following effects: 

(1) Radiation heat transfer from external and internal surfaces of 
the cone 

(2) Convective heat transfer to the cone 

(3) Conductive heat transfer along the cone in the longitudinal 
direction 


In the derivation of K, accomt was taken of the fact that the 
Reynolds number was sufficiently small for the flow to be in the slip-flow 
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station was the same for five different Mach mmibers tested between 
4.4l and 4.59, and observation shows the axial temperature gradient 
along the cone to be negligible. The trend in the recovery factor 
(figs. 6(b) and T) was found to be similar to vhat established in other 
cold-flow tests (see ref. 2) j that is, the recovery factor is essentially 
independent of Mach number and Reynolds nuniber, and for laminar flow can 
be predicted by the sq^uare root of the Prandtl number . Both the low 
value of the local Reynolds number, Re < 2.9 X11 ^ , and the close 
agreement of the experimental data and the •/Pr rule for the recovery 
factor in laminar flow. Indicate that the flow over the cone was 
definitely laminar. 

A 


Hot -Flow Tests 


The results of the hot-flow tests are presented in flgirres 8 
through 11. In figure 8, distributions of surface temperatvire along 
the model are shown for each test condition. The time history of the 
surface temperature is given in fig\ire 9 for cne representative test 
condition. The recovery factors calculated from the basic data of 
figures 8 and 9 are presented in figures 10 and 11 as functions of 
stagnation temperature and Reynolds number (based on T' ) , respectively. 

The recovery factors shown in fig\rres 10 and 11 are in fair agree- 
ment in magnitude with the sqimre root of the Prandtl number. However, 

because of the spread of the d ata, no unique agreement with either 

or n/pt^ can be established ('/Pi^ “ O.85) . The data presented in 
figure 10 also indicate two consistent trends^ one of increasing recovery 
factor with decreasing pressui^e^ and one of decreasing recovery factor 
with increasing ter^jerature. These trends will he discussed in the 
following paragraphs* 


Effect of pressure on recovery factor *- Ihe trend of increasing 
recovery factor with decreasing pressure X^'ta^p.at ion teirperature constant) 
seems to he in agreement with the present thecries* The over-all effect 
of lowering the Reynolds nnmher from the range associated with continuum 
flow, through the transition and slip-flow re^;imes and into the free- 
molecule -flow regime is to raise the recovery factor gradually from its 
continuum value (r ^ 0,85^ assuming laminar f.:.ow) to its free— molecule- 
flow value (r 1*0, see ref. 8) . For the pr^dsent investigation, a 
reduction of stagnation pressure at constant stagnation ten^erature 
lowered the Reynolds number* Since the Reyno.i.ds numhers are in the 
slip-flow range (see fig. , sn Increase in T.he recovei*y factor with 
(j^gQPeaslng pressure would he expected. Also, it was ascertained that 
the assimiptlons involved in the data reduction have little effect on 
this trend. 


00 H 



no H 00 


13 


Effect of teag^eratiure on recovery factor .- The other trend apparent 
in figure 10 is that the recovery factor decreased with increasing 
temperattire. This trend is contrary to what might he expected in the 
slip-flow regime, as discussed previously, since the Reynolds number 
decreased with increasing temperature (constant pressure) , From the 
data of the present investigation, however, it is difficult to ascertain 
the tme cause of this trend. The data may represent the actual effect 
of tenperature on recovery factor at higher temperatures or an apparent 
effect resulting from uncertainties in the experimental measurements 
or the assumptions involved in reducing the data. To see the effect of 
the various factors that may contribute to this result it is necessary 
to review briefly the basic nature of the Investigation and to consider 
the assumptions involved in the data reduction. 

The determination of the recovery factor (given by eq. (l)) is ■ 
dependent upon a knowledge of the wall enthalpy under conditions of 
continuum flow and no heat transfer. For investigations at elevated 
temperatures wherein the model surface becomes heated anti thus radiates 
energy, it is necessary either to cancel the radiation or to formulate 
an energy balance equation to determine the adiabatic wall enthalpy hp. 

Since the latter approach was taken, the value of the actual 
convective heat transfer was determined (all the terms except 

dQctp in the energy balance equation (B21) could be evaluated from the 

experimental data) . The classical continuum convective heat transfer, 
then, was obtained from 


dQc = — ^ (8) 

vhere | is a variable and accounts for various effects of the low<« 
density high— teiirperature flov field. In particular^ for this investigation^ 
t included the effects of changes in boundary- layer thickness^ velocity 
profile^ Maxwell reflection coefficient and thermal accommodation 
coefficient. 


Once the classical heat transfer is known^ the recovery factor in 
terms of the actual heat transfer^ is expressed as 


r 





c 


( 9 ) 


where dQ^^^/^k* is equal to K (see eq. (2)) and is given by equation 

(B2i+) . (The term k* is a function of the geometry and the flow field . 
properties, and is, in effect, related to the classical heat-transfer 
coefficient.) 
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Since the assumptions involved in evaluating the term | can affect 
the value of the recovery factor, it is worthwhile to see what effect 
reasonable variations in these assumptions have t5)on the recovery-factor 
data presented herein. The first assumption involved the boundary- layer 
thickness. The boundary- layer thickness was evaluated by an empirical 
eq^uation, (BI9) . If this thickness were thinner by either I5 or 75 
percent from the value used herein (a reducticn in thickness could be 
produced by either slip-flow effects or the use of a different base for 
the Reynolds number) , the recovery at the high tengjeratures would be 
decreased by approximately 1 and 5 percent, respectively. 

Secondly, the velocity profile of the boundary layer was assumed 
to be quadratic. Changing the order of the profile from second to first 
or to third produced virtually no change in the recovery factor for the 
range of Knudsen numbers of the test. 

The Maxwell reflection coefficient, d, and the thermal accommodation 
coefficient, a, were taken as O.9O for these analyses. A variation of 
-10 percent (from O.90) for d and a producet a 1 - and 2-percent 
increase, respectively, in the recovery factoi'. Very little is known, 
however, of the correct value of these paramei ers for the surface condi- 
tions and wall temperature of the present investigation. A discussion 
of these factors is given in reference 8 where they are indicated to be 
a futnction of wall temperature. This is espec-.ially important for the 
thermal accommodation coefficient since, at the wall temperatures herein, 
the internal vibrational energy was excited aid this may have caused a 
much greater variation in d than considered. 

Also, the accuracy of the basic experimeiital data could affect the 
recovery factor presented herein. The accuracy of the thermocouples for 
the wall measurement was taken as ±1 percent, the manufacturer quoted 
accuracy. A 1-percent variation in wall temperature produces a 
corresponding 1 -percent variation in recovery factor. The stagnation 
temperature was taken as that indicated by the probe. If the probe 
itself had a recovery factor and radiation losses, this value would be 
low, or if the pebble bed radiated to the probe, this value would be high. 
A il-percent variation in stagnation temperature produced a +l-percent 
variation in recovery factor. 

In the preceding paragraphs the factors affecting the behavior of 
the recovery factor in the high-ten^jerature low-density range have been 
discussed briefly. Although no explanation cfin be given for the cause 
of the trend noted herein, it is evident that no single source or combin- 
ation of sources of error of the magnitudes m'jntloned was able to alter 
appreciably the trend of decreasing recovery factor with increasing 
temperature. However, it must be remembered that the correct values of 
some of the parameters are unknown, especially the thermal accommodation 
coefficient, and sufficiently low value co\ild con5)letely alter the present 
trend. Corroborating these results with others must be postponed until 
other high-temperature data become available. 
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CONCLUDING REMARKS 


A 


Wind-tunnel tests were conducted to determine whether the recovery 
factor can he predicted hy the simple square root of the Prandtl number 
^e at stagnation temperatures approaching those encountered in flight. 

number of at a stagnation temperature 
of about 6 to R and for a Mach number of 6.0 at stagnation temperatures 
between 12TO and 2600 R. The Reynolds nuniber for all tests was less 
than 2 . 5 x 10 and for the Mach number 6.0 tests in particular, the 
Refolds number was always less than 10^, thus introducing slip-flow 

GII 0G"US * 


From the data shown, the Prandtl number rule is believed to be 
valid for estimating the recovery factor in the range of temperatures 
considered. The Prandtl number used herein is based on a reference 
temperature T' ; however, because of the nature of the experimental data, 
no conclusion can be drawn as to the necessity of using this tengjerature , 
as opposed to say the static temperature of the gas at the cone surface. 


, recovery- fact or data of the hot-flow tests also Indicated two 

definite trends: one of Increasing recovery factor with decreasing 

pressure and one of decreasing recovery factor with Increasing temper- 
ature. The first trend was seen to be in agreement with the present 
theories; that is, in slip flow an increase in the recovery factor with 
decreasing pressure (or Reynolds number) TOuld be expected. The second 
trend, that of decreasing recovery factor with Increasing stagnation 
te^erature, was shown to be contrary to what might be expected in the 
slip-flov regime since the Reynolds number decreased with increasing 
temperature. The cause of this trend could not be explained, but, it 
was pointed out that no single source or combination of sources of 
error of the magnitudes considered was able to alter appreciably the 
l^end of decreasing recovery factor with increasing temperature. 

trend should be viewed with reservation until other 
high-temperature data become available for corroboration. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Mar. 20 , I 96 I 
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APPENDIX A 


A THEORETICAL SOLUTION OF INVISCID CONICAL FIDW OF REAL 
GASES AT TEMPERATURES BELOW DISSOCIATION 


A necessary prerequisite to the detemilnatton of the transfer of 
energy from a flowing gas to a Body in contact with that gas is a 
thorough knowledge of the properties of the f.-OW field. For flows of 
air wherein the local static temperatures are low^ the air can he 
considered ideal, and convenient solutions exist for most simple flows. 
However, for local static temperatures greater than about TOO R, air 
can no longer he considered an ideal gas and it Is necessary to develop 
an exact solution for the conical flow of a real gas. 

Because the solution presented herein is of interest to the general 
field of high-temperature gas dynamics, it is presented in a manner to 
maEe it entirely independent of the experimental portion of this paper. 

The solution, as discussed in this part, has been programmed for 
an IBM 70^ computing machine and is available upon request to Ames 
Research Center. 

Determination of the properties of a conical flow field involves 
the simultaneous solution of the following equations J 


(pv sin 0) + 2pu sin 0 = 0 

(Continuity) 

(Al) 

(30 



du 

(30 

( Irrctationality) 

(A2) 

udu + vdv+^=0 
P 

(Momentum) 

(A3) 

h + "7^^ ' 

(Energy) 

(At) 

h = h(p,p) 

(state) 

(A5) 


A suitable combination of equations (Al) , (AP) and (A3) yields the 
well-known Taylor-Macoll equation for conical flow (see ref. l8) 



- ^ 


A 

3 

1 

8 




d^u 

de^ 


(A6) 



rO H GO 


2D 


17 


which is "basic to the solution presented herein. The use of the above 
set of equations assumes that the gas is inviscid, is everywhere in 
thermodynamic equilibri-um, and is isentropic before and after the shock 
wave . 


The eq\iation of state for the gas used in this paper is the Beattie- 
Bridgeman equation with provision for variable specific heats. Choice 
of this model for the gas was based on the following: 

(1) This equation of state is an accurate representation of gas 
characteristics up to temperatures vhere dissociation (diatomic 
and polyatomic molecules) and ionization (monatomic molecules) 
are encountered. 

(2) This equation is accurate at densities where intermole cular 
forces should he considered. 

(3) Constants required in the Beattie-Bridgeman equation of state 
are available for a large nxmiher of gases. 

Table I presents the constants needed to make the Beattie^-Bridgeman 
equations applicable for a number of the more common gases of current 
interest. For a more complete table of the Beattie-Bridgeman constants 
see reference I9. Data pertaining to the vibrational energy levels of 
other gases can be obtained in references 20 and 21. 

The discussion of the solution is divided into three main areas; 

(1) Determination of flow characteristics ahead of the shock wave. 

(2) Determination of flow characteristics immediately behind the 
shock wave. 

(3) Determination of the flow field between the shock wave and the 
body. 


DISCUSSION 

Conditions Ahead of Shock Wave 


In general^ the conditions ahead of a shock wave are defined by one 
of two sets of three variables. In free-f light Investigations the 
variables usually specified are velocity^ static temperature, and static 
pressure. In wind-tunnel investigations the variables usually specified 
are Jyiach number, stagnation temperature, and stagnation pressure. The 
first set is used in this paper as the direct input to the calculation 
of the flow through shock waves. Hence, where the second set is specified, 
a considerable amount of calculation is required to transform these 
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variables into those of the first set with the inclusion of real-gas 
effects. The required transformation equations are the svibject of this 
section, and as presented herein are solved t y an iterative process. 

In solving for the temperature and pressure that satisfy the set 
of initial conditions, recourse is made to the usual conditions of 
adiabatic reversible flow. However, to accoint for the differences 
between an ideal gas and the Beattie-BridgeniEn gas (see ref. 22), the 
enthalpy, the entropic equation of state, and speed of soimd equations 
are modified as follows: 

h = Jh = <IcpT (A7) 


7 


Pt Pt Kt W w 

(A8) 

a = n/i" a = \/ I7RT 

(A9) 


The correction factor for the enthalpy eqixation is given by: 


J = 


1 


+ 


7-1 

7 


G + p 



2Ao ^ 
RT " ’jgJ 


+ 0^ ^ 
V2 RT 2 


Bob 




+ P" 


\ T^ /J 


where G Is given below In general form and the necessary constants 
are presented in table I for the gases considered of interest herein. 


G 



^i^i^vp 


exp(0vi/T)-l- 


i = 1.2,3> 


The correction factor for the entropic eqviatlon of state is given 

by: 


K = 


E 

F 


(l + eip + C2P^ + 63p3) 


(AlO) 


where e, E, and F are given below (see table I for necessary constants) 


no H CO 
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G 


Aof 


“ M' ■^O 


- B^b - 


BnC 


©3 = 


BqBg 


E = 


ex^D 




evi/T 


.exp(0vjyT)-lJ 
1 - exp(-0vi/T) 


diXi 


1 = 1,2,3,^ 


F = exp 




The correction factor for the speed of sound is given hy: 


I = 2 (l + 2eip + 3eaP^ + ^eap^) 

7 


(All) 


In order to calculate the ahove correction factors^ the density 
and specific heats of the gas must be calculated for the Beattie-Bridgeman 
gas. The density is given by 

P ^ ^ (l + glP + g2P^ + gsP^) (A12) 


vhere 


gl = - ^ 


ga = 


o 2 

2ex ^es 

(rt)^ 


[ -Qs] 

(rt)^ 


gs “ 
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The specific heat at constant volume is given hy: 


6Rc 


Cy Cy + “^vih iji3 P 1 2 




Bnh 


2 


(A13) 


where given helow and Cy as well e.s the constants for 

are given in table I. 


^vih ^ ^ ^ %^i 
i 


^Vi/ 2T 


sin h(0Yi/2T), 


The specific heat at constant pressure is given hy 




= c^ + R 


1 + 


||p)(l+BcP +Bohp2^ 


1 + 2exP + 362P" + ^esp3 


The ratio of specific heats^ y , is hy definition 

• n 


(Al4) 


(A15) 


Since the forms are known for the varlouj: correction factors and 
thermodynamic properties for the Beattle-Bridgeman gas^ the iterative 
process to determine the static temperature^ static pressure^ and velocity 
proceeds in the following manner. As a first approximation^ the static 
temperature and pressure are taken as the ideal-gas values; namely^ 


A 


T = Tt 



(Al6) 


P = Pt 




(A17) 


The values for static ten^ierature and pretssure are then used to 
calculate the correction factors for the Beattie-Bridgeman gas. The 
static temperature for the Beattie-Bridgeman gas is then found from 

T = (l-X) (A18) 
J 


po H CO 
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where 


X = 




tjf + £ 

7-1 I 


The static pressure is fotind from 


P = Pt 



_2L 

7-1 

Kt 


(A19) 


and the velocity is given hy 


U = Ma 


(A20) 


The degree of accuracy of the ten^jerature and pressure given hy 
equations (A 18 ) and (AI 9 ) is dependent upon the initially assumed 
values given hy equations (AI 6 ) and (AI 7 ) • It is usually necessary to 
Iterate the solution which Involves using the values calculated in 
equations (AI 8 ) and (AI 9 ) to re-calculate the correction factors and 
then repeat the solution of equations (AI 8 ) and (AI 9 ) . 

The significant characteristics ahead of the oblique shock which 
are required to calculate the flow immediately behind the oblique shock 
are the vectoral quantities 

vectoral quantities perhaps bears some explanation. When the momentum 
and energy equations normal to the wave are written 

P„ + p„U ^ = const = (A21) 

00 '-'00 


hr 




con 


= const = h+ = h(T+ , P+ ) 
^“n “n ^ 


(A22) 


it is observed that the normal stagnation conditions must be used. These 
quantities in general are different from the free-stream Mach number and 
stagnation conditions. The normal Mach nimiber is given by: 

= ^toSln 8s (A23) 

The evaluation of the stagnation temperature and pressure normal 
to the shock "wave proceeds as follovs^ where the static conditions and 
Mach number are given in eq^uations (A 18 ) ^ (AI 9 ) and (A23) ^ respectively. 
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The stagnation temperature normal to the shock wave is given hy 

TdoJoo 


Tt 

W J. 


toon^l-X) 


(A24) 


where X is evaluated, for The stagnatjon pressure normal to the 

shock wave is given hy 


PKtco. 




n 


(A25) 


K 




A 


As a finsi approximation to the ahove temperatures and pressures^ 
the ideal gas values of 


and 


V) 

= P» + ^ M»o^) 7-1 


(A26) 


(A27) 


can he used where T and P are the free-stream values determined pre- 
viously ( (A18 ) , (AI9) ) . This solution is also Iterated until sufficient 
agreement is found. 


Conditions Immediately Behind Siock Wave 


If the conditions in front of the shock wave are known, the static 
temperature and pressure behind the shock wave can be found by an 
iterative solution. As a first approximation, the following values 
can be used: 


Tg = Ik, 


-(7-1)] [(7-1)1^,^^+ 2 ] 
(7+1) 


and 


Ps = Pco 


r 2 mo/ -(7-1)1 


(A28) 


7+1 


(A29) 
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A 


The density hehlnd the shock wave is given hy the following: 



(A3O) 


where 


L - 


+ giP + 


The ten^De rat lire and pressure for the Beattle-Bridgeman gas are 

then: 

Tg = ^ (l-Xs) (A3I) 

^ s 

where 

Xs - ^^32) 


and 


PsLs(^’"^s) 


(A33) 


Here agaln^ the required agreement for successive values of Tg 
and pg Is obtained hy it e rat ion - 

The total Mach number and velocity behind the shock wave are given 
by the following relation. 

The normal Mach number behind the shock is: 


The tangential Mach namiber is given hy 

Ms, tan = ^,te:nj ^ ^^35) 



2k 


and the total Mach number is then given by: 

Ms = A/ + Ms^,tan (A 36 ) 

The speed of souQd hehind the shock vave is given hy: 


ag = s/ IfRTs 


and the total velocity is 


U — Mgag 


(A37) 

(A38) 


The stagnation conditions behind the shock wave are found from: 


Tfc 


s 


^s'^s 

Jts(l-Xs) 


(A39) 


A 


where X is given by equation (A32) and from 


Ptg 


Ks VTtgy 


(Ako) 


The ideal-gas flow equations are used to determine the initial values 
of Ttg and ptg used in the iterative solution of equations (A39) 

(A 40 ) . 


Conditions Between Shock Wave and Body 


In the conical flow fields between shock wave and body^ the velocity 
and hence temperat-ure and pressure are constant along radial vectors 
from the cone apex. The solution of the conical flow field is then a 
step-by-step solution from one radius vector to the next. 

The velocity along any radius vector is given by; 


X = ^i-i 


d0 


de + 


i-i 


-) 

d0^. 


d0' 


1-1 


2 


(AM) 
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The velocity normal to the radius vector is given by: 





(a42) 


A 

3 

1 

8 


where the Taylor-Macoll equation is used as follows; 



(ug + cot e 




(A6) 


The subscript i denotes any given radius vector^ the first vector 
(at an angle 0g-Z^) being i=l and so forth. The shock wave would 
then be i=0^ at which the following boundary conditions are known: 


0 = 0s 

u = U^cos 0s (a4-3) 


The entire solution is completed by numerical integration at each 
subsequent radial vector (0-A0) until the radius vector equals the cone 
surface half angle. The total velocity on any given radius vector is 
given by 


U = s/u2 + v2 (A44) 

At this point in the solution^ values for stagnation temperatiire 
(A39 ) } stagnation pressure (a4o) , and velocity (Akk) have been obtained. 
Hence ^ with the exception of Mach number ^ the same quantities are known 
as were used in determining the flow characteristics ahead of the shock 
wave. As a first approximation to the value of Mach number^ the following 
is used: 


Me = 


Ue 


a,Q 


i-l 


(a45) 


This value of Mach number and the stagnation temperature and pressure is 
then used in equations (A18) and (AI 9 ) to calculate the static tempera- 
ture and pressure. Again this solution must be solved by iteration 
because of the approximation for Mach number. 
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A question that always arises is^ what quantity^ if any^ can he 
used as a convenient and invariant normalizing parameter for the velocity. 
For this puipose^ the only quantity that is truly invariant is the 
stagnation enthalpy^ which is analogous to the square of the stagnation 
speed of sound of ideal- gas theory^ and is a parameter used to normalize 
the velocity 


W = — ^ {Ak6) 

The solution is continued until Vq = 0^ which is the boundary 
condition for the cone surface. A 

In using a solution of this type^ for any assumed values of 
and the half-angle of the cone is unknown, and given only by the 

complete solution. If the conditions for a cone of definite half -angle 
are desired^ several solutions can be run^ and the flow field for the 
cone interpolated between two solutions. Shown in table II are sample 
solutions for 10^^ 20°^ and (half -angle) cones at 1000^^ 2600^^ and 
3600^ R. 

The largest effect of increased temperature on the properties of 
a conical flow field is seen to be in the values affected by the internal 
energy of the gas (i.e.^ vibrational energy) . The enthalpy is naturally 
expected to Increase with temperature^ and the calculations show^ for 
the higher temperatures^ significant changes in the specific heat^ Cp^ 
and 7. The pressure and^ hence^ the pressure drag are seen to be 
unaffected by temperature. 

For the tenperature and cone angle considered^ the effects of 
increased temperature are seen to be small. However^ for larger cone 
angles, or higher tenperatures, considerable departure from the ideal-gas 
conical flow properties exist. 


no H CO 
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APPENDIX B 

DEVELOPMENT OF EQUATION FOR RECOVERY FACTOR 


As pointed out in the introduction^ several authors have found 
that for flight speeds greater than a Mach number of 3*5 the enthalpy 
rather than the temperature of the gas must he used in evaluating the 
recovery factor. For such conditions then^ the recovery factor is 
defined as 

^r*"^c 



The quantities h^ and h-^ vere determined directly from the measure- 
ments of the flov characteristics of the wind tunnel^ and the properties 
of the conical flov field vere determined hy the methods discussed in 
appendix A. The quantity h^ is the enthalpy of the air at the surface 
of the cone under the conditions of no heat transfer to or from the cone 
and of te335)eratures at equlllhrium. The first condition did not exist 
during the present test. Furthermore, the skin temperature of the cone 
rather than the enthalpy term^ h^-, vas measured. Hence ^ it is necessary 
to use an involved set of equations to reduce the skin-temperature data 
to the desired enthalpy quantity, h^. The development and discussion of 
these equations is the subject of this appendix. 

To develop the eq-uations for the recovery factor, it vas necessary 
to determine the heat balance into and out of an incremental volimie of 
the cone in the vicinity of the ten 5 )erature -measuring station. The 
incremental volume, dV, shovn in the sketch belov is bounded by tvo 
sides having areas, A^ and A^, and the external and internal siirfaces 
of area, Ag and A^^ respectively: 





T 


X 
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To the order of the differential lengthy dx^ these areas and volnines 
can he expressed in terms of as follows 


As - Ai + ^ dx 
dx 

(B2) 

^i “ (^1 “ Jtm^) — — 

^ m cos rp 

(B3) 

dAg = (Ai + rtm2) — — 

^ m cos p 

(Bk) 

dV = A^dx 

(B5) 

Ax = Jt(2xm tan cp - m^) 


dA. 

-T— = 2mii tan cd 
dx 



The various heat flows comprising the net total heat to the volume^ dV^ 
are: 

1. The conductive heat transfer along the metal skin^ Qi and Qq* 

2. The radiative heat transfer from the cone surfaces^ and Qre* 

3* The convective heat transfer from tie air to the cone surface^ Q,^. 

Each of these forms of heat transfer will he discussed in the following 
sections . 

It should he noted that during the present tests ^ the cone tempera- 
tures used to determine the recovery factor were measured after the 
model reached equilibrium temperatures so that steady- state techniques 
were used in the following derivations* 


CONDUCTIVE HEAT TRANS! ER 


In deriving the conductive heat-transfer equations, two assumptions 
were made: 

1* No temperature gradient across the skin thickness. 

2. No circumferential temperature gradient. 


OD H Oo > 
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The equation for the conductive heat transfer into dV is given hy: 

dT 

Qi = ^ (B6) 

and the heat transfer out is given by 

Qq - (B7) 

^ (Note: The minus sign arises in equations (b 6 ) and (B7) because x is 

in a direction of positive heat flow; that is^ Ty decreases with positive 
so that dTy/dx will always be negative. The minus sign must 
therefore be included for these equations to be compatible.) 

The net transfer of heat by conduction to volume dV is then given by 
(neglecting second-order terms in dx) 


^^oi “ ^i “ Qo 


= !<■ — 


dA d.Ty 
dx dx 


+ An 


d^T 


w 


dx' 


dx 


(b8) 


where A2 has been expressed in terms of A^; km is the thermal 
conductivity of the metal; and Tw is ten5)erature measured at the wall. 

The above equation requires a knowledge of the temperature gradient 
along the cone. In determining the values for dTy/dx and d^Ty/dx^^ the 
data for each thermocouple station were fitted to a second-order 
expansion of the form T ~ ax^ + bx + c. In this way a consistent 
approach to obtaining the derivatives was assured. 

Also required in the above equation is the thermal conductivity of 
the material. For the material used (stainless steel A16) the conduc- 
tivity was obtained from reference 23. The variation with temperatiore 
of the conductivity is shown in figure 12. 


RADIATION HEAT TRANSFER 
Radiation to External Surroundings 


For the present investigation^ the radiation heat transfer was 
determined experimentally by allowing the model to cool in the wind 
tunnel under conditions of zero flow, a near vacuum, and tunnel-wall 
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temperatures the same as those used to measure the recovery factor. The 
ten^erature of the cone as it cooled was measured as a function of time 
and the q_uantity (d'IV/dT)zF relation to the cone temperature vas 
obtained and is shorn, in figure 13 • During these tests to determine 
the radiation heat transfer^ it tos noted that the teii^eratures at the 
three measuring stations along the cone surface and at the front 
radiation shield vere the same^ thus indicating no conductive heat 
transfer along the cone or internal radiation to the shield. Hence ^ 
the radiation heat transfer from dV to the external surroundings vas 
given by: 


dQre 





Substituting equation (B 5 ) into the above equation gives 


A 


<iQre 





(B9) 


where 

Pjjj density of material 

Cpjjj specific heat of material 

T time 

Unfortunately, experimental data giving the specific heat of the 
material, Cpjjj, as a function of temperature for stainless steel hl6 
could not be found in the literature. Consequently, the values used 
for Cpjjj were derived from the empirical expression: 


°Pm ~ 



Ci 

Mi 


(BIO) 


where 

Wi fraction of specie 1 
Q± molar heat capacity of specie 1 

M;j^ molecTlLar weight of specie 1 
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The chemical con 5 )Osltion of the stainless steel is given in table III 
and values of cp^^^ in relation to temperature are shovn in figure l4. 

(As a matter of interest other values of and Cp^ for different 

steels are presented in figures 12 and l4, respectively. Data are 
from references 2^ to 28.) Note that the experimental Cp values for 

stainless ^03 (coniposition nearly the same as kl6y see table III) has 
a similar trend to that given by the theory for ^l6 and ^ 03 . 


Radiation to Internal Shields 


In contrast to those measurements made to determine the external 
radiation vhile the model vas cooling^ the temperature of the cone 
surface and the internal radiation shields vere not the same for the 
steady state (flov-on) tests. The equation for thermal radiation 
between the internal surface of the cone and the radiation shield was 
derived using the approach of reference 2^ and is given by 


dQrj_ - ^d(T^ - Tg"^) (Fi-s + Fs*-i)dAi 


where 

G emlssivlty of unoxidlzed platinum 

cf Stefan-Boltzmann constant 

Tg temperature of front radiation shield^ *^R 

Evaluating the form factors gave: 

^i-s == sin 9 
^^s-i = 0 

and if equation (B3) is substituted for dAp^ the internal radiation is 
given by: 

- Gd(V - Ts4)sln qp dx (B11) 

The emlssivlty;, e, is shown as a function of temperature in figure 
15* The data were obtained from reference 29* 
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CONVECTIVE HEAT TRANSFER 


In this section, the basic equation for convective heat transfer 
to a flat plate in continuum flow will he presented. Various correction 
terms to the basic equation will then be presented in order that the 
results be applicable to the present investigation. 


Basic Equation 


The equation for convective heat transfer to a flat plate incorpo- 
rating the variation of the fluid properties with temperature is given 
■by ( see ref . 4 ) 

dQc = 3600 CHPcUc(br ~ l^r)dAe 


where 

Cg dimensionless heat-transfer coefficient 
hr recovery enthalpy 

h^ enthalpy at wall 

air density at edge of boundary layer 
Uc velocity at edge of boundary layer 

and the factor 36OO converts the unit time from seconds to hours. 

It will be noticed that the convective heat transfer is then pro- 
portional to the difference between the recovery enthalpy, which is 
required to determine the recovery factor, anl the enthalpy at the wall, 
which can be con^^uted from the temperature measurements made in the 
inve st 1 gat i on . 

Van Driest (ref. 4 ) calculated the dimensionless parameter, 

Cg 7 RSq? ^ laminar boundary layer on an 1 isulated flat plate at 

various stagnation teii^jeratures . His result Ls plotted in figure I6. 

Expressing Ag in terms of Ai, the basic eq-uation for the convec- 
tive heat transfer becomes 


, V Pc^c / V -h 

dQe = 3600 (Cg^) (hr - hw) ^ ^os cp ' 

V Ren 


(B12) 
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Correction Factors Considered 


The basic equation is for an insulated flat plate in continuum 
flov vhereas in the present case^ a cone in the slip-flow regime is 
being considered (see fig. 4) . Correction factors are req_uired to 
account for these differences. A summary of the theoretical and 
experimental approach of the effect of slip flow on heat transfer is 
given in reference 30 the applicable equation is 


A 



(B13) 


where the bar denotes slip flow and B is a function of the slip 

velocity ratio. In addition^ the transverse curvature of the cone will 

have an effect on the skin friction. In reference 30 or 31 it was 

shown that when 5*cos cp/R is less than or of the order of unity (for 

the present tests 6*cos cp/R^ where 5* is boundary- layer displacement 

thickness and R is local radius of external surface of cone^ was 

always less than 0,30) , the skin- friction coefficient^ cp ^ is given by: 

xr 

Cftr = 
or 

Gftr = 

where the subscript tr denotes transverse curvature. 

Combining the correction factors and considering the effect of 
transverse curvature gives the convective heat transfer in the slip-flow 
region: 


dQ 


ctr 






(B15) 


where the factor transforms the basic equation from that for a 

flat plate to that for a cone (see ref. 3^) • 

The factor B is given by (see ref. 30) : 


B == 




- c, 


1 _ 


1 - C2 


% 

u, 


(Bl6) 



3k 

where 


87 f 1 - ^ cr \ 

^ " (7 + 1) \ «■ J - V 


C2 = 1 




\ 


7 + V 




n Kn 

— = — — • (slip velocity ratio) 

and 

n order of velocity profile 

c_ 

7 ratio of specific heats^, — t- 

Cv 

cj Maxwell reflection coefficient 

a thermal accommodation coefficient 

Kn Knudsen number 

Pr’ Prandtl number evalmted at reference temperature 

The factor G is given by (see ref. 3I) : 


T, 


G = 1 + I 


0.517 + 0.913 ^ + 0 . 121 ( 7 - 1 )Mc" 


A 


1 / 

tan (p 


vhere 
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Ce 


l^c 


and 

1^ coefficient of viscosity with silbscripts w referring to wall and 
c to outer edge of "boundary layer 

The coefficient of viscosity, p, is given hy Sutherland's equation 
(ref. 12 ) as 


A 

3 

1 

8 


m3/ 2 

li = 264.084x10“^ gQ - g - 

The factor ("cf/cf) is given hy (see ref- 30) 

^ 5-5[f(L)]^^ 

(2 + L) 

vhere 



(BIT) 


(B18) 


and 


f(L) 




1.0 


16 - 

15(2 + L) 


The Knudsen number is required to evaluate several of the previous 
terms. The Knudsen number as used herein is defined as: 



vhere the local mean free path^ as determined from kinetic theory 
is given by 
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and the homdary- layer thickness for an insiilated flat plate (ref. 33) 
is given hy 

5 =. (1 + 0.08 y ^) (B19) 

By the application of the Hantsche -Wendt (ref. 3^) transformation 
from a flat plate to a cone the Knudsen nimiber can he expressed as^ 


Kr 


0.756 Me 
(1 + 0.08 Mc^) 



(B20) 
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Correction Factors Neglected 


The effect of honndary- layer dlsplacemert on heat transfer vas 
neglected Because it wa_s found that for each test condition the hypersonic 
interaction parameter [X = Cq/s] Re q) ] never exceeded a value of 0.3 

(see ref. 30) . Such lov values iniply that the Interaction Between the 
Boundary layer and Invlscld flow region is weak. 


ENEBGY BALANCE 


Since the data were taken at equlliBriuE. conditions^ the net heat 
transfer to the volinne dV can Be written as 

0 - dQoi - dQre - dQrp + dQc^.^ (B21) 


^The transformation gives for the cone values 

^c - ^/T (see ref. 3^) 

5c = B//3 (see ref. 35) 


where A and 5 are the flat-plate values. 
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Siibstltution of egmtlons (b 8) , (B9) , (B11) and (BI5) gives 


0 = km ^ dx + PmCpmAi 


d^Tw' 


dx dx 


dx^ 


dT, 


V 


dr 


dx 


ZF 


-eo(T„‘ - Osin cp * '111 




3600 


r I \ Pc^c /, , \ (At + Jtm^) 

(GHn/R^) (hr “ hvj-) ^ ^ 


sf Bee 


dx 

(B22) 


Thus, it Is possible to write 

= IV + K (B23) 

where K is the total correction to the measured enthalpy of the air 
at the wall to determine the recovery enthalpy. 

In terns of the measured quantities , 



where Cp is the specific heat at constant pressure. 
Solving equation (B22) for K gives 


K = 


6236 BG 


Cf 


P(,U^(Ai + Jtm^) 




-1 


-k. 


dAj 
dx dx 


kuiA 


d^T. 


V 


^ dx^ 




T.^)tan cp (^" 


m 


(B2A 


in which all terms can he calculated from the measured cone temperatures, 
wind-tunnel flow characteristics, and the relationships derived herein. 
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TABLE I.- GAS CONSTANTS FOR BEATTIE-BRIDGEMAN EQUATIONS 
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TABLE III.- COMPOSITION OF STAINLESS STEEL 4 l 6 AND 403 


Element 

Composition, 

percent. 

416 

403 

Carloon^ C 

0.15 

0,15 

Chromium Cr 

13.00 

13.00 

Iron^ Fe 

83.97 

34.85 

Manganese, Ma 

1,00 

1.00 

Molybdenm, Mb 

0.60 


Nickel, Ni 

— 

0,50 

Phosphorus, P 

0.04 

— 

Silicon, Si 

1.00 

0.50 

Sulfur, S 

0,24 
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Figure 2.- Model and instrumentation details 
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Figure 5-“ Transport and thermodynamic properties of air as a function of teirperature. 
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(c) Specific heat at constant pressure 
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Figure 5»- Continued. 
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Figure 5.- Concluded. 
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Figure 6.- The temperature ratio and recovery factor as a function of the distance along the model. 
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Figiore 6.- Concluded. 
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Figure 7,- Recovery factor as a function of local Reynolds n-umber for the cold-flow model. 
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Figure 8.- Temperature distribution along hot-flow model. 
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Figure Representative surface ten^^erature-time history for hot-flow 

model . 













Thermal conductivity, k 



Figure 12.- Thermal conductivity of 

a function of 
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iveral types of stainless steel as 
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Figure 13 •- Rate of change of surface temperature with time as a function 
of surface temperatiu-e for the radiation tare runs (Mqo = 0, = O) . 
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Figure l6.- Local heat-transfer parameter for a laminar boundary layer on an insulated flat plate 

as a function of stagnation temperature. 




